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We present an analytical derivation of the Sachs Wolfe effect sourced by a primordial magnetic 
^- , field, generated by a causal process, such as a first order phase transition in the early universe. 

As for the topological defects case, we apply the general relativistic junction conditions to 
match the perturbation variables before and after the phase transition, in such a way that the 
total energy momentum tensor is conserved across the transition. We find that the relevant 
contribution to the magnetic Sachs Wolfe effect comes from the metric perturbations at next- 
to-leading order in the large scale limit. The leading order term is strongly suppressed due 
to the presence of free-streaming neutrinos. We derive the neutrino compensation effect and 
confirm that the magnetic Sachs Wolfe spectrum from a causal magnetic field behaves as 
1(1 + 1) Cf oc I 2 as found in the latest numerical analyses. 
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The origin of the large scale magnetic fields observed in galaxies and clusters is still unknown: 
one of the possible explanations is that they have been generated in the primordial universe. A 
primordial magnetic field of the order of the nano Gauss c ould leave a detectable imprint in the 
cosmic microwave background (CMB) anisotropics^ ' ' ' '. Here we concentrate on its effect on 
the temperature CMB spectrum at large scales and more particularly on the Sachs Wolfe effect. 
The motivation is that conflicting results exist in the literature regarding the ^-dependence of 
the Sachs Wolfe effect induced by a causal primordial magnetic field: the analytical analysis 
ofW found £{£ + 1) Cf scaling as l~ l or more negative, and the same result was found in the 
numerical calculation of^ on the other hand, I ' ' ' found £{£ + 1) Cf scaling as £ 2 . The aim of 
this paper is to explain this discrepancy analytically (for a more detailed derivation see ). 

We assume that a magnetic field is generated in the early universe by a sudden phase transi- 
tion, as for example the electroweak (EW) phase transition. We consider a stochastic magnetic 
field with no background component and we suppose that the magnetic energy momentum tensor 
is first order in perturbation theory. We study the effect of the magnetic field on the metric and 
fluid (matter plus radiation) perturbations by solving analytically Einstein's and conservation 
equations in the long wavelength limit. We take into account the neutrinos in our derivation. 
In order to connect the solutions before and after the magnetic field generation, we match the 
geometry and the fluid variables at the phase transition time, so that the induced three metric 
and the extrinsic curvature are continuous . This implies the conservation of the total energy 
momentum tensor across the phase transition, and it completely determines the metric and 
fluid perturbation variables after the magnetic field generation. Before neutrino decoupling, we 
find that at leading order in the large scale expansion k/7i <C 1, the metric perturbation $ is 
proportional to <3? oc H-q{J-L / k) 2 , where IIb is the magnetic field anisotropic stress. This ind uces 
a contribution in the CMB spectrum scaling as £(£ + 1) Cf ex £~ l , and consistent with^^. 
However, once neutrinos decouple and start free-streaming, they acq uire a non-zero anisotropic 
stress, which acts to compensate and reduce the magnetic field one^^. We demonstrate that 
this compensation drastically reduces the leading order contribution to the CMB spectrum, and 
that the dominant contribution becomes the one from the n ext-to - leading order in the k/H <C 1 
expansion, which induces then £(£ + 1) Cf oc £ 2 , as found in ^ l 4 I ^ I. 



2 Solutions for the metric and fluid variables 

In this section we solve for the metric and fluid perturbations. We consider only scalar perturba- 
tions on a spatially flat Friedmann background, and we work with gauge invariant variables. We 
consider a stochastic primordial magnetic field with spectral index n > 2 because of its causal 
generation^: (Sj(k)S*(q)) = (27r) 3 5(k — q)(<% — kikj) Ak n , where A is the amplitude of the 
spectrum. We work under the one-fluid MHD approximation, meaning that the conductivity of 
the universe is high, so that we can neglect the electric field. The magnetic field is characterised 
by its energy density p-g, anisotropic stress 7Tb and Lorentz force is that describes the momen- 
tum exchange between the magnetic field and the primordial fluid. These quantities satisfy the 
relation ps/2 = vtb + 3^b/(2&), coming from momentum conservation. To solve for the metric 
and fluid perturbations, we combine the system of Einstein's and conservation equations into a 
second order differential equation for the gauge invariant variable D, corresponding to the total 
(matter plus radiation) density perturbation in the velocity-orthogonal slicing. The equation is 
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where a dot denotes derivative with respect to conformal time rj, np is the anisotropic stress of 
the matter plus radiation fluid, Ob = Pb/ Pr&A an d IIb = ^n/Praxi- To determine the Sachs Wolfe 
effect, we need to solve this equation for scales which are over the horizon at recombination. 
Therefore, we can drop the terms proportional to (kfH) 2 . 

We split the universe's evolution in three stages: before magnetic field generation, between 
magnetic field generation and neutrino decoupling, and after neutrino decoupling. Before mag- 
netic field generation, the universe is filled with matter and radiation only, for which we con- 
sider adiabatic initial conditions. Since at early time neutrinos are coupled to photons, their 
anisotropic stress is zero. Consequently, the two last lines of Eq. (jT]) vanish and we solve the 
homogeneous equation for L>LLUJ. At time tjb , the magnetic field is generated by a causal process 
that we assume to act 'fast', i.e. within one Hubble time, as for example a sudden phase transi- 
tion. After tjb, the second line in Eq. ([1]) acts as a source for D. We solve it using the Wronskian 
method. In order to connect the solutions before and after the field generation, we match the 
variables in such a way that the total energy momentum tensor is conserved across the transi- 
tion. Following what has been done in^ for the analogous case of the topological defects, we 
find that this matching imposes the continuity of the variables <j) and V at the phase transition 
(see Eqs. (36) of^). Finally, at time rj v , neutrinos decouple and start to free stream. They 
acquire a non-zero anisotropic stress *k v that acts as a new source term in Eq. ([!]): 7Tf = Rv^u 
(where R u = p u /Pra,d = 0.4). Solving Eq. ([1]) using the Wronskian method requires to know the 
time evolution of the neutrino anisotropic stress, ir u . The evolution of n^ in the presence of an 
external constant anisotropic stress has been studied in^ (see also^). The neutrino anisotropic 
stress quickly adjusts to the external one and compensates it (see appendix B or^for a detailed 
derivation of tt u ). The final time dependence of tt^ is rather complicated (Eq. (B.12) of ) , but 
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is related to the amplitude of the primordial potential. With this approximation, we can solve 
analytically for D. 



3 Sachs Wolfe effect at leading order in the long wavelength expansion 

With our solution for D valid through the whole universe's history, we compute the metric 
variables $ and ty and the velocity perturbation V, using Einstein's equations. We then split 
the fluid into its individual components and solve for the photon density perturbation D g7 and 
velocity perturbation V 7 , using standard adiabatic initial conditions. With this we compute the 
Sachs Wolfe contribution to the temperature anisotropy at large scales, given by (see e.gMty 
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where x Tec = kr] TCC <C 1 for superhorizon scales at recombination. This contribution is pro- 
portional to the magnetic field generation time (e.g. ye — 1CP 12 for generation at the EW 
phase transition), and to the neutrino decoupling time, y v ~ 10 -6 . It is therefore strongly sup- 
pressed: f{y r cc-,yv,y'B) — 10~ 6 . This contribution corresponds to the effect of the magnetic field 
anisotropic stress from its time of generation to the neutrino decoupling time. The subsequent 
magnetic contribution to the Sachs Wolfe effect, arising from y u up to recombination time, is 
cancelled by the free-streaming neutrinos. With the above result we compute the CMB spec- 
trum. Neglecting the integrated Sachs- Wolfe, and using that, for a causal magnetic field with 
n > 2, the power spec trum of the magnetic field anisotropic stress is constant in k up to the 
damping scale ko ' ' ' 

(n B (k)n B (q))=(27r) 3 5(k-q)|n B (A0| 2 with |n B (£0| 2 =n^^-^, wefindEl (3) 
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where fffaec) den otes the visibility function. Therefore, we find that the C MB spectrum scales 
as l/£, as in^^, rather than as I 2 as found in' 1 ^'. The reason why refs.^™ do not find the 
l/£ dependence is because the initial conditions that they insert into their Boltzmann code are 
derived after neutrino decoupling, when the magnetic field anisotropic stress has already been 
compensated by the neutrino one. Eq. @ shows that the period of time between magnetic field 
generation and neutrino decoupling leaves an imprint on the CMB spectrum. This imprint is 
however much too small to be observable. Therefore, we now proceed to compute analytically 
the next order contribution in the kfH <C 1 expansion, that will lead to a £ 2 -dependence in the 
CMB spectrum. 

4 Next-to-leading order contribution to the Sachs Wolfe 

In order to compute the next-to-leading order contribution to the temperature anisotropy, we 
need to solve for the metric variables at next-to- leading order in the k/% <C 1 expansion. The 
easiest way to compute this order is to use the curvature perturbation £ = ~ < ^ + 3 ( 1 2 _ lu \ y& — y\- 
Indeed, starting from the leading order solution that we computed for D, we can calculate the 
next-to-leading order in the curvature Q. We use then this solution to compute the next-to- 
leading order in <I>, by integrating from the definition of £. This can be done analytically if 



"Note that the magnetic field does not affect the curvature at leading order, 0(l/x^ ac ). It is therefore still 
conserved on large scales. 



we approximate the time evolution of the neutrino anisotropic stress at next-to-leading order 
as ir v (y) = (dif^e + ^IleXy — y v ) 2 x 2 ec , where d\ and d,2 are two arbitrary constants that we 
determine from the conservation equations of the neutrino and photon fluids. With this, we find 
for the temperature anisotropy at next-to-leading order 

— (Mo)- -0.2 ft B -2.7IIb • (5) 

The energy density spectrum and the cross-correlation one have been calculated irJ ' ' \ and 
they share the same /c-dependence as the anisotropic stress spectrum, Eq. ([3]). Denoting Ct the 
amplitude of the energy density spectrum and C the cross-correlation one, we find for the CMB 
spectrum 

£(£+l)Cf ~ g 2 (?? r ec) [o.04fi + 7.29 5 + 0.54 c] ^J- i( f + ,\ ■ (6) 
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Therefore, we confirm that the next-to-leading order contribution to the Sachs Wolfe effect scales 



as £{£ + 1) Cf oc £ 2 , as found numerically inl^l 

5 Conclusion 

In this work we present an analytical computation of the Sachs Wolfe effect induced by a 
primordial magnetic field. We have restricted our analysis to a magnetic field generated by a 
causal process, such as a first order phase transition. In this case, the initial conditions for 
the metric and fluid variables are determined unambiguously by imposing conservation of the 
total energy momentum tensor across the transition. Using these initial conditions, we have 
computed analytically the leading order and next-to-leading order magnetic contribution to the 
Sachs Wolfe effect. We have found that the leading order contribution is sourced only by the 
magnetic field anisotropic stress, and leads to a CMB spectrum scaling as l/£. However, this 
contribution is strongly suppressed once the magnetic field anisotropic stress is compensated by 
the one of the neutrinos. As a consequence, the dominant contribution to the Sachs Wolfe is the 
next-to-leading order one, that generates a CMB spectrum scaling as £ 2 . Our analytical work 
solves therefore the discrepancy regarding the ^-dependence of the magnetic Sachs Wolfe. 
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